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Abstract 

We evaluate the probability that a loop of string that has sponta- 
neously nucleated during inflation will form a black hole upon collapse, 
after the end of inflation. We then use the observational bounds on the 
density of primordial black holes to put constraints on the parameters 
of the model. Other constraints from the distortions of the microwave 
background and emission of gravitational radiation by the loops are 
considered. Also, observational constraints on domain wall nucleation 
and monopole pair production during inflation are briefly discussed. 



1 Introduction 



It has been proposed |], 0] that loops of string can spontaneously nucleate 
during an inflationary period of expansion in the early universe. This is a 
quantum tunneling process, somewhat analogous to the spontaneous nucle- 
alion of spherical bubbles of true vacuum in the problem of false vacuum 
decay 0. It has been shown that loops can nucleate at considerable rates 
provided that their tension fi is not much larger than if 2 (here H is the 
expansion rate during inflation) and, in this case, the distribution of loops 
emerging from inflation may have significant cosmological consequences. 

In this scenario, loops nucleate with sizes of the order of H" 1 , and are 
subsequently stretched to large sizes by the inflationary expansion. After 
inflation, in the radiation dominated era, the loops eventually fall within the 
Hubble radius and start oscillating under their tension. 

Because the instanton describing the nucleation of a loop is maximally 
symmetric, the nucleated loops tend to be nearly circular. It is well known 
that an exactly circular loop would form a black hole upon collapse [f|, §]. 
However, the nucleated loops will not be exactly circular. The reason is that, 
during inflation, quantum fluctuations on the string on subhorizon scales gen- 
erate shape instabilities as their wavelengths become larger than the horizon 
[§, |7| (see also ||). The question then arises of how many of the nucleated 
loops will be circular enough to form a black hole after the end of inflation. 
This is the problem we address in the present paper. 

The plan of the paper is the following. In Section 2 we calculate the 
probability of black hole formation from a nucleating string. In Section 3 we 
use this result and the observational bounds on the abundance of primordial 
black holes to put constraints on the parameters of the model (which are 
essentially fi and H.) Additional constraints from the distortions in the 
microwave background and the emission of gravitational waves by the strings 
are considered. In Section 4, we briefly discuss observational constraints on 
domain wall nucleation and monopole-antimonopole pair production, which 
can also occur during inflation |l] . We summarize our conclusions in Section 
5. Some details are left to the appendices. Appendix A deals with linearized 
perturbations on a circular loop of string in flat space, and in Appendix B 
we consider the effect of damping processes on the probability of black hole 
formation. 
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2 Probability of black hole formation 



In this section we estimate the probability that a loop of string that has 
nucleated during inflation will form a black hole once it collapses under its 
tension, after the end of inflation. 

The classical worldsheet of a circular loop of string that has nucleated 
during inflation is given by |I| 



where R is the physical radius of the loop, H is the expansion rate of the 
inflationary universe and t is the usual cosmological time in the flat Friedman- 
Robert son- Walker coordinates. The parameter to can be interpreted (in a 
loose sense) as the time at which the loop nucleates. Throughout this paper, 
we shall work in the approximation in which the strings are infinitely thin, 
although nucleation of strings whose thickness is not much smaller than H^ 1 
is also possible 0. 

The loops nucleate with a size of the order of the horizon size R ~ H^ 1 , 
and afterwards they are stretched by the inflationary expansion. For R >> 
H- 1 they grow proportionally to the scale factor 



After inflation, in the radiation dominated era, the loops continue to be 
stretched by the expansion of the universe, R(t) oc t 1 / 2 , until they enter the 
cosmological horizon, at t ~ R. Soon after the loop comes within the hori- 
zon, the effects of expansion become negligible and the loop starts behaving 
approximately as it would behave in flat space. 

Let R c be the radius of a particular loop at horizon crossing. The mass of 
this loop is 2tiR c h, where \x is the string tension. The Schwarzschild radius 
corresponding to this mass is 



As the loop shrinks under its tension, its rest mass is converted into kinetic 
energy, so that the total energy of the loop remains constant. It is clear that 
if the loops were exactly circular, they would all eventually shrink to a size 
smaller than R s , thus forming black holes (of course this argument assumes 




(1) 



R(t) » H- l e H{t ~ to) . 



R s = AttG/jR, 



(2) 



3 



that R s is larger than the thickness of the string, which will always be true 
for sufficiently large loops.) 

However, loops nucleated during inflation will not be exactly circular. 
Quantum fluctuations around the symmetric solution (|l|) will cause small 
departures from the circular shape, and this will determine the probability 
of black hole formation. Consider a circular loop of radius R(t) lying in the 
z = plane and centered at the origin. Using cylindrical coordinates (p, B, z), 
we can parametrize a perturbed (non-circular) loop as follows: 

p(0,t) = R(t)+A r (0,t), (3) 

z(0,t) = A t (0,t). 

That is, we decompose the perturbations into a radial part A r parallel to the 
plane of the unperturbed loop, and a transverse part At perpendicular to that 
plane. In eq.@, p and z are physical coordinates (rather than co-moving.) 

To study the evolution of a loop, we shall divide its history into three 
epochs: 

a- The loop nucleates and then expands during inflation. 

b- After inflation, while the loop is larger than the cosmological horizon it 
is stretched by the expansion until it crosses the horizon. 

c- Once the loop comes within the horizon, it shrinks under its tension. 

We will follow the evolution of the perturbations A r , At through the epochs 
a, b and c. For this it will be useful to expand them in Fourier modes: 



A,= E 



L=2 



a( i) cos LB (2) sin LB 

^r,L 7= + ^r,L r- 
\/7T a/7T 



(4) 



and similarly for A t . Note that the sum does not include the modes with 
L = and L = 1. This is because, as it was shown in Ref. [0, such modes 
do not correspond to true perturbations, but only to spatial rotations and 
spacetime translations of the unperturbed solution ([!]). 

a- Evolution during inflation. In the process of nucleation, A£ L (i = 1,2; 
A = t, r) have to be treated as quantum variables. The linearized theory of 
quantum fluctuations around the classical solution (0) was studied in Ref. [0], 
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using a covariant formalism. The worldsheet ([I]) has the internal geometry 
of a 1+1 dimensional de Sitter space 0, and it was shown that the pertur- 
bations A r and A t behave as two uncoupled scalar fields of tachyonic mass 
—2H 2 'living' in this lower dimensional de Sitter space. The symme- 



m 



tries of the problem suggest that when the string nucleates, the fields A r , A t 
should be in a de Sitter invariant quantum state (see also Ref. 0). This state 
was constructed in Ref.fTj using the Heisenberg picture. The corresponding 
Schrodinger picture wave functional #({A? L }) can be obtained using well 
known manipulations (see e.g. ||). Since we are dealing with a linearized 
theory, the probability distribution associated with \l/ has the Gaussian form 



&}) 



1^1 



11(2* 



a 



\,L) 



-1/2 



cxp 



\,i,L 



2 a 



where the standard deviations <7a,l are given by |7j 



alAt) =<*|(A«) 2 |*> 



2fiL 



H 2 R 2 (t) 
1? - 1 



+ 1 



(5) 



(6) 



For wavelengths much larger than the horizon, HR » 1, we have o~l oc 
R, i.e., the amplitude of the perturbations is conformally stretched by the 
expansion (recall that R grows like the scale factor during inflation). For 
wavelengths much smaller than the horizon the second term in (|6|) dominates 
and we have o~l ~ constant, as expected. This second term is just the flat 
space contribution to the quantum fluctuations, and we shall subtract it 
in what follows in order to avoid the well known ultraviolet divergences in 
<Ai>. 

Soon after a given wavelenth becomes larger than the horizon, it can be 
treated as classical. Here we adopt the point of view that Eq.fl^) gives the 
probability distribution for the initial amplitudes of the perturbations, which 
will evolve classically thereafter. 

b- After inflation: while the loop is larger than the cosmological horizon. 
As we mentioned before, the loops are conformally stretched until they enter 
the horizon. Since we do not have analytical solutions to describe this epoch, 
we shall use the following approximations |T0| , 11]: for wavelengths larger 



than the horizon the amplitude of the perturbations is conformally stretched 
by the expansion, whereas for wavelengths within the horizon the amplitude 
of the perturbations remains constant. 
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Perturbations with wave number L enter the horizon at time tL ~ -R(^l) / L. 
At this time we have a L = [2fiL(L 2 - l)]-^ 2 HR(t L ). The loop enters the 
horizon at the time 



t c ~ R(t c ) 



R 2 



(tL) 



LR(t L ). 



At this time we still have the same <jl than at tL, since the amplitude of the 
perturbations is frozen after ti, so we write 



a L (t c ) = [2/xL(L 2 - \)]- x f 2 HR[t L ) = [ 2 ^L\L 2 - l^HR,. (7) 



c- Loops inside the horizon. Once a loop comes within the horizon we can 
ignore the effects of the expansion of the Universe and consider the collapse 
of a loop in flat space. In this case the evolution of the perturbations can 
be solved analytically. The details can be found in Appendix A. We find 
that during this period the behaviour of radial perturbations A r is different 
from that of transverse perturbations A t . The amplitude of radial perturba- 
tions shrinks by a factor of L as the loop collapses, while the amplitude of 
transverse perturbations remains constant. For R << R c we find 



Introducing (H) in (|5[), we obtain the probability distributions for the vari- 
ables A^ L at the time when the radius of the unperturbed loop enters the 
Schwarzschild radius R s . 

The dynamics of gravitational collapse and black hole formation from a 
given loop configuration is, of course, a complicated issue that cannot be 
treated in detail analytically. However, we know that if the amplitude of the 
perturbations is much smaller than R s at the time t s when the unperturbed 
loop enters the Schwarzschild radius, then a black hole will form. On the 
contrary, if the perturbations are large compared with R s the loop will prob- 
ably self intersect and fragment into smaller loops before it shrinks to a size 
R s (the resulting loops will be far from circular, so a black hole will never 
form in this case). To compute the probability of black hole formation, we 
shall count a particular loop as a black hole with gaussian weight 



0~ r ,L 



(2/iL 7 ) 




(8) 



l 2 ER c . 




(9) 
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Here 5 is the r.m.s. perturbation when one averages over the circunference 
of the loop (at time t s ) 

Z7T Jo Ztc rr* T 

where we have used In Eq. (0), a is an unknown parameter of order 
1, which models our lack of knowledge on the dynamical details. We should 
emphasize that the precise form of the "Gaussian window" ([5p is not very 
important here; we have taken the Gaussian shape just for computational 
convenience. Had we taken, for instance, a step function window, our final 
conclusions would be the same. 

From @ and @, the probability of black hole formation from a loop that 
has nucleated during inflation is given by 

v bh = j i*i V-*"/**.)- n ^ L = n n fi + 1 . (io) 

J \,i,L \=r,tL=2 \ ' iUC -"-s/ 

Using (H), we can rewrite it as 

oo / T5\- 1 ( /-5\ _1 

= II (! + 2^1 [ 1 + 2 lj) ' ( n ) 

where we have introduced the notation 

L* = [167r 2 5(aG/i) 2 ]- 1/5 . 

Here B = 4nfi/H 2 is the Euclidean action of the instanton describing the 
nucleation of the loop. Note that Vbh does not depend on R c , the radius of the 
loop at horizon crossing. Note, also, that the dependence on the parameters 
if, /i and a is only through the combination L*. 

The physical meaning of L* is the following. Modes with L » have 
o"l << aR s . These modes will always have amplitudes much smaller than the 
Schwarzschild radius and they will not contribute to Vbh- This is also clear 
from ([Til). Therefore is basically the number of modes whose amplitude 
can be large enough to prevent the formation of a black hole. 

The function Vbh{L*) cannot be given in closed form. However, one can 
compute it numerically to arbitrary precision by including a sufficient number 
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of terms in the product (|TT|). A plot of Vbh versus is given in Fig. 1. From 
the graph we see that the probability decays (essentially) as an exponential 
function of L*. 

In deriving Eq.(|ll]) we have neglected damping processes, such as gravi- 
tational radiation and friction. In principle these could reduce the amplitude 
of the perturbations and therefore increase the probability of black hole for- 
mation. In Appendix B we briefly discuss the damping mechanisms and the 
limits in which they can be neglected. We find that damping can be safely 
ignored from the calculation of Vbh provided that Gji is in the range 

{a 2 B) -2/9 < Gfl<5 . 1 -V5) 1/3 , (12) 

where m p is the Planck mass and M is the mass of the black hole. The upper 
limit is due to gravitational radiation, which is important for heavy strings, 
while the lower limit corresponds to friction due to the surrounding matter, 
which is important for light strings. In the next Section we will be interested 
in Vbh for black holes of mass M ~ 10 19 m p , so Eq. flT2D leaves us with a wide 
range of values of Gfi for which ( ITT] ) is valid. 



3 Observational constraints on loop nucleation 

During inflation, loops of size ~ H" 1 are produced at a constant rate per unit 
physical volume. After nucleation these loops grow like the scale factor, so we 
expect that the number of loops with size ~ R contained in a volume ~ R? 
will be independent of R. That is, we expect a scale invariant distribution 
of loops . The number density distribution of loops will be given by []T[ 

dN dR 

w = v w (13) 

(Actually, the distribution ( |T3| ) has a lower cut-off at R ~ H" 1 , since we do 
not have loops smaller than that. It also has an upper cut-off at R ~ EH~ L , 
where E is the total e-folding factor since the onset of inflation. See Ref. [|IJ 
for details.) 

The coefficient v is the number of loops produced during an expansion 
time if -1 in a volume H~ 3 . This has been estimated in the semiclassical 
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approximation in Ref. JI|, using the instanton methods || 



v = Ac' 3 . (14) 

Here, as in Section ^ B = Att^/H 2 is the euclidean action of the instanton 
describing the nucleation of the string (which is just a spherical worldsheet 
of radius H~ l ). The coefficient A has not been calculated and we shall leave 
it as a free parameter here, presumably of order one. Let us now consider 
various observational constraints that one can place on the parameters of the 
model. 

a- Constraint from primordial black hole abundance. 

After inflation, the distribution of loops with R » t is still given by (|T3|), 
since these loops continue to be stretched by the expansion. Loops with size 
smaller than the cosmological horizon, R « t, are not stretched and are 
simply diluted by the expansion, 



dN _ 

dV ~ v 



a(R) 
a{t) 



where a is the cosmological scale factor. For loops that enter the horizon 
during the radiation dominated era, we have 

dN (l^dR 

W =V \Tr) It' (16) 

Multiplying by the mass of the loops, M pa 2 r n[iR 1 we obtain the mass density 
distribution 

/2ttG7^ 3/2 



w-teJ dM - (17) 

Dividing by the critical density p c = 3/327rGt 2 we have, at the present time, 



V ; p c 3 V w \GMJ M ' V ; 

where t eg ~ 4 • 10 10 sec. is the time of equal matter and radiation densities. 

As we described in Section 0, some of the loops in the network will form 
black holes with probability given by eq.(|ll|). To obtain the spectrum of 
primordial black holes that are produced through this mechanism we multipy 
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fll8|) by Vbh- The fraction of the density parameter Q in black holes of mass 
~ M is then given by 

n*(M) « n,^5(^) 3/ v (i^) V2 . (19) 

The strongest observational constraint on the abundance of primordial 
black holes comes from the emission of 7-rays by holes that are evaporating 
at the present time [TJ|]. These black holes have a mass M ~ 5 ■ 10 14 g, and 
the constraint is given by 

Q bh (5 ■ 10 14 2)£l(T 8 . (20) 

From (|1^) we have 

P bh ■ (G>) 3/2 ^10- 28 . (21) 
We can rewrite this inequality as 

- log 10 V bh (L*) - ^ log 10 (aG/x) + (log 10 e)B^N, (22) 

where N = 28 + log 10 (v4a:~ 3 / 2 ). Note that the left hand side only depends 
on B and the combination aGfi. Therefore, for given N, eq. (|22]) will exclude 
a certain region in the plane (B,aGfi). Since the parameters a and A are 
expected to be of order 1, their contribution to iV will be small. Note also 
that even if the observational constraint (p0|) was improved by one order 
of magnitude, this would only increase N by one. At any rate, we expect 
iV m 28 plus or minus a few units. 

In Fig. 2, the boundary of the excluded region in parameter space is 
depicted for iV = 27, 28 and 30. It is seen that the result is not very sensitive 
to the value of iV and, quite generically, the bound will be satisfied provided 
that 

aG>~10- 4 , (23) 

independently of B. This can be easily understood since, for such values of 
aGfi, the parameter L* will be very large and, correspondingly, the proba- 
bility of black hole formation will be exponentially small (see Fig. 1). 

Also, from Fig. 2 we see that for large values of B, say -B~50, the 
constraint on G\x is practically removed due to the exponential suppression 
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in the number of nucleated loops. Of course this limit is not very interesting 
cosmologically. 

b- Constraints from the microwave background 

A string moving across the sky produces a discontinuity in the observed 
temperature of the Cosmic Microwave Background (CMB) between both 

roughly of order (ST/T) « 10G/i. For values of 
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sides of the string 

v ~ 1, the loop sizes extend all the way to the present Hubble radius, and 



the bound obtained in Ref. 14 



G//~10" 



[v ~ r 



should apply. However, for v « 1 this bound should not necessarily be sat- 
isfied, since in this case all the loop radii are much smaller than the horizon. 

Let us first consider the effect of loops with R > t eq . The angle 9 sub- 
tended by the string is given by 



y = tan(6/2) 



R 

d 



1 - 



dp 
3t 



(24) 



where t is the age of the Universe and d is the present distance to the string. 
Loops of given radius start oscillating when R ~ t, and therefore they will 
only have an effect on the CMB provided that they are at a sufficiently low 
redshift. Using the fact that the light that is now reaching us from a point 
at a distance do was emitted at time t = t [l — (do/3t )} 3 , loops that were 
oscillating when the observed light passed through them must satisfy 



R<t 



dp 
3t 



(25) 



The number of loops with radius in the interval [R, R+ AR] at a distance 
in the interval [do, Ad Q ] is given by 



AN(R, d ) 



v 



RH\ 



AirdlARAdo- 



(26) 



Changing variables {R, d } — > {y, d } through eq. ( p4[ ) and integrating over 
all possible distances d allowed by the constraint (|25f) we obtain 



dN(y) = 36nu 



1 




dy. 



(27) 
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The expectation value of the number of loops subtending an angle larger 



than = 2 tan 1 y is then 
N >& = / dN[y] = 







In 






1-3 


K 


hi 


) 


y 











This function is plotted in Fig. 3. It is seen that for > 10° (the angular 
resolution of the COBE experiment) and v < 10 -3 , for instance, the expected 
number of events iV>e is only of order 1. This means that for such low values 
of v we can hardly expect to see any loops that large. 

However, one must take into account that 'unresolved' loops whose an- 
gular size is smaller than the detector beam angular size can still produce a 
signal, although of lower magnitude. For simplicity, in what follows we shall 
assume that the loops are nearly circular. In that case, if the plane of the 
loop is perpendicular to the line of sight, the temperature distortion in the 
CMB has the profile of a 'top hat' with radius equal to the radius of the 
loop at the time when the light rays crossed this plane ||13|| . On the angular 
scale set by R (the radius of the loop when it is at rest), the temperature 
fluctuation is \5T/T\ ~ 10G/i (the sign depending on whether the loop is 
expanding or contracting at the time when the light rays cross the loop). On 
larger angular scales, the signal will be inversely proportional to the solid 
angle that we are considering. 

Therefore if the detector beam has an angular size j3 larger than the 
angular size of the loop 0(-R, do), the measured temperature fluctuation will 
be 

1 5T /0\ 2 

A "io^T S! U)- < 28) 

For small angles, the distribution ( |27| ) can be approximated by dN(Q) = 
487rz/0 _3 (i0, and using (|28|) to express as a function of A we have 

247TZ/ 



iV >A 



Here N > & is the expected number of events due to unresolved loops that 
give a signal larger than 5T/T on a detector with beam angular size (3. This 
equation can be rewritten as 



uGji 



240tt T ' 



(29) 
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where (3 has to be expressed in radians. 

The COBE experiment constrains 5T/T to be less than 10~ 5 on scales 
(3 pa 10° pa .17 rad. It is not clear what is the maximum number of events 
allowed by the COBE data. Taking N&~1 seems a little over restrictive 
since one small string can easily be hidden behind the galaxy. To be more 
conservative we shall take N&~10, which from fl29|), yields the bound 



(v < 1(T 2 ) (30) 

This bound comes from the contribution of unresolved loops, and it only 
applies for v < 1CT 2 . For larger values of z/, the number of 'resolved' loops is 
sufficiently large that the usual bound 

G/i-lO" 6 [v > 1(T 2 ) (31) 

should be used. 

So far we have considered loops with R > t eq . A similar analysis can 
be carried out for R < t eq to show that these loops cannot produce relative 
temperature fluctuations in excess of 10~ 5 provided that the constraints (|30), 



fl3~T|), and Eq.(p3|) from the previous subsection are satisfied. 
c- Constraint from gravitational radiation 

The millisecond pulsar observations place a constraint on the density 



parameter in gravitational radiation of period ~ 1 year (T5| , 

VL g <4-l(T 7 /r 2 . (32) 



From strings, we have [16 



1287T (GA l/2 

where Q r ~ 4 • 10~ 5 /i~ 2 is the present density parameter in radiation and 
7 9 ~ 100. With this, the bound (|52|) reads 



KG7i) 1/2 ^2-10- 3 . 

This is always satisfied provided that the bounds (^),(^(]) and (|3~1~D from the 
previous subsections are satisfied. 
c-Structure formation 
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A question of cosmological interest is whether the constraints derived so 
far are compatible with a scenario in which the nucleated strings would seed 
the formation of large-scale structure in the Universe. Let us make a rough 
estimate of the values of the parameters needed for such purpose. 

Assuming, for simplicity, cold dark matter, the mass accreted by a loop 
of radius R < t eq at the present time is 

M = 2irfiRz eq . (R < t e q) (33) 

Here z eq ~ 2 • 10~ 4 h 2 is the redshift at time t eq , with h the Hubble constant 
in units of 100A;m s~ 1 Mpc~ 1 . This is because perturbations start growing 
at time t eq and they grow proportionally to the scale factor. From (|T6| ) the 
number density of loops of radius ~ R is hr ~ u(t eg R)~ 3 ^ 2 z~ q . Using (|33|) , 
the mean separation between objects of mass M is given by 

j -1/3 -1/3 ( Mt eq Z eq \ ^ 

dM~n R( * M) ~v I y-^^J ■ (34) 

Rich clusters of mass M c \ ~ \^h~ x M & are separated by distances of 
order d c i ~ 50h~ 1 Mpc. Therefore, from (|34"D , we obtain the normalization 

z/ 2 / 3 G/i « 5 ■ Kr 8 /^ 1 . (z/^10- 2 ) (35) 

This normalization is valid for G/i~10 _6 /i (hence z/~lCU 2 ). For smaller values 
of G\x the loops that accrete masses of order M c i have R> t eq , and (^) does 
not apply. 

For G//~10 _6 /i, loops start accreting matter when they enter the horizon, 
and we have 

M = 27T^R (JJ . (R>t eq ). 
Using tir = u/(Rto) we obtain 

dM ~ n'll^ ~ z/~ 1/3 Af (27T//)- 1 . 
Again, normalizing for rich clusters, we find 

v^GfM w 2 ■ 1(T 7 . (z/^KT 2 ) (36) 
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The region excluded by the constraints fl23|), fl30|) and ( pTj ) in parameter 
space is plotted schematically in Fig. 4 (shaded region). It is seen that the 
parameters needed for structure formation (thick solid line) [from ( |3"5j ) and 
(36[)1 lie marginally in the allowed region. Of course, the normalizations (|35|) 
and (|36D should not be taken too literally, and one expects large error bars 
in the thick solid line of Fig. 4. However, a scenario in which the nucleated 
strings may seed some of the observed structure in our universe does not 
seem to be ruled out. 

We should now mention the effects of compensation [17]] which has been 
ignored in the above discussion. When a loop nucleates, its energy is bal- 
anced by a corresponding deficit in the local densities of matter and radiation. 
This compensation reduces the gravitational effect of the loop on surround- 
ing matter and on background photons on scales greater than R. When a 
compensated loop radiates away its energy, it leaves behind an underdense 
region. As the density contrast grows, this region will evolve into a void with 
a dense central object seeded by the loop. The compensation is somewhat 
reduced by radiation and neutrinos streaming into the underdense regions 
and can be further reduced by loop fragmentation. The effect of compen- 
sation on structure formation has not been fully investigated, but it seems 
reasonable to assume that this effect will not be dramatic on scales crossing 
the horizon before t eq . Then the relations (|35|),(|36D should still be valid by 
order of magnitude. 



4 Domain walls and monopoles 

Just as in the case of strings, spherical domain walls can spontaneously nu- 
cleate at a constant rate per unit volume during inflation. As a result, the 
Universe will get filled with a scale invariant distribution of walls given by 
( P~5| ) and (|T4D, where now B = 2ix 2 aH^ z is the action of the instanton de- 
scribing the nucleation of the wall. 

The observational constraints on domain wall nucleation are very different 
from the ones that one can impose on strings, since the mass of a spherical 
wall 

M = AnaR 2 (37) 
grows quadratically rather than linearly with the radius. An inmediate con- 
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sequence is that walls of size 



R c > (SnGa)- 1 (38) 

will all collapse to form black holes KJ, since for them R s > R c . 

To estimate the microwave background anisotropy induced by domain 
wall bubbles, one has to take propper account of the compensation effect 
discussed at the end of the previous section. This is a somewhat complicated 
issue and we shall not attempt to address it in detail here. However, one can 
make a rough estimate of what this effect should be by using the following 
arguments 

In a pure dust universe the compensation would be complete, and back- 
ground photons would be unperturbed on scales greater than the co-moving 
scale of R s at horizon crossing. Therefore, on larger scales, the whole effect 
should be due to the underdensity in radiation that is needed to compen- 
sate for the mass of the black hole. This underdensity causes radiation from 
neighboring regions to move in, partially destroying compensation. As this 
happens, the initial underdensity propagates away from the black hole as a 
sound wave at the speed of light. At time t eq the mass of the underdensity, 
M r , is comparable to the mass of the black hole, M b h- At later times, we have 
M r « Mbhz/z eq , due to cosmological redshift. The gravitational potential of 
M r extends up to scales of order t, and we have <p ~ GM r t~ x . This will cause 
temperature distortions of order 

5T J GM bh z 5 / 2 



T to z eq 

on scales of the order of the cosmological horizon at redshift z. 

The angle a subtended by the horizon at redshift z is given by tana w 
(z 1//2 — so at small angles 5T/T oc a -5 . In particular, for a ~ 10°, we 
have 

ST M bh 

— ~ — — (a « 10 ), (39) 
T M u ( J 

whereas for a ~ 90° we have a much lower effect 5T/T ~ M bh /(M u z eq ). Here 
M u = G^ta is approximately the mass of our observable universe. 

Let us denote by R ma x the radius of the largest wall that may have existed 
in our observable universe (up to distances d = 3£ ). How massive can this 
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wall be? From fl39|), and using 5T/T~10 5 , we have 

M max ^10- 5 M U » 10 5 V (40) 

From (|37D and ( fi"0| ) we can see that -R maa; must be less than £ e? even if a 
is as low as the electroweak scale. Using 



dN 



^ dR, 



dV R s/H% 2 V t 
the number of walls with size larger than R within a volume 47r/3(3to) 3 is 



iV >i? = 247TZ/ 



Therefore, the larger wall that we can expect to find in the observable uni- 
verse has radius R ma x ~ [247r^t e g 1 / 2 t ] 2//3 , and Eq. fl4~0|) yields the bound 

^/3^_< 1() -64 (41) 

Let us now turn to monopole pair production during inflation. Unlike 
walls and strings, monopoles are not stretched to enormous sizes by the 
inflationary expansion. They are just diluted. An inmediate consequence is 
that the only monopoles that will be relevant at the end of inflation are the 
ones that have been produced during the last expansion time. This density 
is given by U 

n = H 3 Ae- 2nm/H , (42) 

where m is the mass of the monopole (here, as in the previous sections, the 
exponent 2nmH~ l is the action of the relevant instanton.) 

A constraint on the abundance of nucleated monopoles comes from the 
fact that they should not recreate the monopole problem that inflation was 
aimed to solve. Of course this can be easily achieved by choosing the ratio 
m/H to be sufficiently large, but as we shall see, one does not necessarily 
have to impose that. 

Actually, the number density of monopoles at the time of reheating de- 
pends very much on the details of how inflation ended. Consider, as an 
example, a model in which the usual exponential expansion is followed by a 
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short period of power law expansion a <xt p that starts at time t\ and ends at 
time t 2 {p may be larger that 1, in which case we have power law inflation, 
but it need not be). At time t 2 reheating is completed and we enter the usual 
radiation dominated era. 

It is clear that, if no monopoles are created during the power law epoch, 
we will have, at time t 2 , 

nit,) ■ /* (43) 



where n(t\) is given by 



n(t 2 ) 
and 

HI 



T 2 

1.6Af 1/2 - r 



Here M is the effective number of massless degrees of freedom at reheating 
and we have used t\ = pH^ 1 and t 2 = 2p m p (10A/')~ 1 / 2 T~ 2 , with m p the 
Planck mass and T r the reheating temperature. Even if we assume that 
monopoles continue to be produced at a rate per unit volume given by 

AH e ~, 

where 7i = a /a is the instantaneous expansion rate, it is not difficult to 
show (for 2imiH~ 1 > 3p — 4) that the density of monopoles during the power 
law era is dominated by the ones that were already present at t%. Therefore 
Eq . (|43|) is still valid in this case. 

After reheating, the monopoles continue to be diluted as a -3 whereas 
the density in radiation decays as a -4 . Therefore the contribution of the 
monopoles to the density parameter at the time t eq ~ 10 10 s given by 



a, 



m 



n(t 2 ) /t eg \ 1/2 



Pc{t 2 ) 



10 29 /3[P - 



H 



1 3/2 



Ae 



27T171 111 

H 

m„ 



The Parker bound |T8| requires fl 
for the monopoles), so we must impose 

3/2 



~10" 



(m/m p ) (assuming Dirac charges 



H 



L m PJ 



f\p-2)Ae 



2irm <-" 

~10" 



33 



As mentioned above, this bound can be trivially satisfied by supressing the 
creation of monopoles with a sufficiently large m/H ratio. Alternatively, if 
p > 1/2, then Q m will decrease with the reheating temperature through the 
supression factor /. For instance, taking p = 3.5 and / ~ 10~ 4 the bound is 
automatically satisfied without any assumptions on m/H. 



18 



5 Conclusions 



In this paper we calculated the probability of black hole formation by a string 
loop spontaneously nucleated during inflation. The calculation was based on 
the quantum theory of perturbations on strings in de Sitter space developed 
in Ref.[[5J. The result is given in Eq . (|TT|) and in Fig. 1. We then used the 
probability ([11]) to derive observational constraints on the mass parameter 
of the string and on the loop density parameter v. The strongest constraint 
comes from the emission of 7-rays by evaporating black holes. Additional 
constraints are due to the absence of characteristic hot and cold spots that 
would be produced by large oscillating loops on the microwave sky. The 
gravitational radiation background generated by the loops imposes no new 
constraints. The excluded region of the parameter space is sketched in Fig. 
4. 

We have also briefly discussed the cosmological implications of nucleated 
domain wall bubbles and monopole-antimonopole pairs. For sufficiently large 
bubbles, the formation of black holes is inevitable. Although these black holes 
can have enormous masses, their effect on structure formation and on the 
background radiation anisotropy is diminished due to the compensation effect 
[O. The predicted density of the monopoles is given by Eq.([|). Depending 
on the values of the parameters, it can be too high, negligible, or it can lie 
in the interesting range and be potentially observable. 
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Appendix A 

In this Appendix we study linear perturbations to a circular loop of string 
in flat space. Before going to the particular case of a circular loop, we shall 
consider perturbations to an arbitrary string configuration whose worldsheet 
is given by x M (^ a ) (£ a are arbitrary coordinates on the worldsheet). This will 
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be a trivial extension of the general results given in Refs. || |7]] for the case 
of domain walls. 

Following Ref. M, we parametrize the perturbed worldsheet x M as 



2 

= x <*(£) + Y, n A y A , (Al) 

A=l 

where n Afi are the two vectors normal to the unperturbed worldsheet. They 
satisfy 

n^daXp = 0, n A »n* = S AB . (A2) 

The perturbation fields (p A have the meaning of normal displacements to the 
worldsheet, as measured by an observer that is moving with the unperturbed 
string. 

The dynamics of the string is governed by the Nambu action 

S[x»(£)) = -nfy/=g d% (A3) 

where g is the determinant of the metric induced on the worldsheet g a b = 
daX^dbX^. The equations of motion that result from (|A3|) are well known 

Ox^iC) = 0, (A4) 

where □ stands for the covariant d'Alembertian on the worldsheet. Multi- 

i A 



plying (|A4j) by n A and integrating by parts one obtains 



9 ab K A = 0, (A5) 

where K A b = —d a n A dbX^ is the extrinsic curvature corresponding to the 
normal n A . 

The effective action for the perturbation fields <p A on a given background 
can be obtained by introducing ( |A1|) into the action ( |A3| ) and then 
expanding to quadratic order in (p a . After some lengthy algebra, the result 
can be written as 

Six"] = S[x»] + S . 

The first term is just the action for the unperturbed solution, while the 
second is given by 

1 



d% (A6) 
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with 

In deriving this effective action we have used the equations of motion ( |A5|) 
to eliminate the terms linear in <p A . From (|A6|) we see that, in general, the 
two fields <p A are coupled to each other in a complicated way. 

Particularizing to the case of a circular loop, matters will simplify con- 
siderably. The unperturbed worldsheet in cartesian coordinates (t, x, y, z) is 
given by 

x^ = (t, R c cos 6 cos — , R c sva.9 cos — , 0). (A7) 
R c R c 

Here R c is just the radius of the loop at t = (when the loop is at rest). The 
two vectors normal to the worldsheet, can be chosen as 

= ?-( srn tt ) cos 9, sin 9, 0), 

^ cos R c 

nf = (0,0,0,1). (A8) 
The first one is a radial normal vector, while the second one is transverse to 

(2) 

the plane of the loop. It is easy to see that in this case K ab = and 5 = 0, 
so the fields (j)^ and (f)^ decouple from each other. They simply behave like 
free scalar fields in the curved geometry of the unperturbed worldsheet. 

In particular, the transverse perturbation behaves like a massless mini- 
mally coupled field. The corresponding equation of motion 

□0( 2 ) = -L=d a {V=gg ab d b <P {2) ) = o (A9) 

can be solved trivially. The metric on the worldsheet is given by 

dsl = g ab di a di b = cos 2 {t/R c )[-dt 2 + R 2 c d9 2 }. 
Introducing this metric in (|A9|) we obtain a flat-space wave equation 



(2) - ^0" (2) = 0, (A10) 

where a prime denotes derivative with respect to 9. Eq. ([A10| ) is not sur- 
prising, since in 1+1 dimensions the conformal coupling is the same as the 
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minimal coupling (see e.g. ]nj). The mode solutions of ( |A10| ) are just 
standing- waves that oscillate with constant amplitude. Note that (j)^ 2 ' co- 
incides with A r of Section ^], so we have the following result: transverse 
perturbations to a circular loop oscillate with constant amplitude as the loop 
collapses. 

The equation of motion for the radial perturbations (f>^ can be found 
from (|A|) (with K$ =S = 0). We have 



□0W + #£ ) #< 1 ><V 1 > = o. 

The extrinsic curvature can be obtained from (Wft) and (|A8|) : 



(All) 



U E> ' 

K 



K, 



(i) 



(i) 
to 



0. 



The proper perturbation <ffo\ as measured by a local observer that is 
moving with the string, is related to the radial perturbation A r defined in 
Section through a Lorentz contraction factor 



l \l-R 2 . 



Writing ([All ) in terms of A r and decomposing in Fourier modes [as in 
we have, after some algebra, 

x 2 / t \ • (L 2 - 1) A 



The general solution to this equation can be written as 

Al(0) 



L 



R 2 Lt R Lt 
1 — — sm — + L— cos 



R r 



Rr 



+ 



Al(0)- 7 '' 



L 2 



\ 



1 



r 2 



cos ■ 



R r 



T R . Lt 
L— sm 



R r 



(A12) 



where Al(0) and Al(0) are the values of the perturbation and its derivative 
at t — (we have dropped the subindex r), and R = R c cos(t/R c ) is the 
radius of the unperturbed loop. 
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In the cosmological problem that we are interested in, t = corresponds 
to the moment when the loop enters the horizon and starts collapsing. The 
probability distribution for the initial conditions A^(0) will be given by 
and (0). We also need the initial conditions A^(0). When the wavelength of 
a perturbation is within the horizon we expect < A| >= (L/R c ) 2 < A| >, 
where the brakets indicate average over one oscillation period. Therefore, by 
the time the loop enters the horizon, the perturbations will have developed 
velocities of order A(t = 0) ~ (L/i? c )A^(t = 0). Using this in ( |A12| ) we have 



A L (R « R c ) « Al{Rc 



L 

That is to say, the radial perturbations shrink by a factor of L as the loop 
collapses. 



Appendix B 

Damping mechanisms, such as gravitational radiation and friction, may de- 
crease the amplitude of the perturbations on a circular loop of string and 
therefore increase the probability of black hole formation. In this appendix 
we briefly discuss the limits in which the effects of gravitational radiation 
and friction can be ignored. 

Gravitational radiation can only smooth out perturbations on wavelengths 
smaller than X g = ^gGfit, where 7 9 ~ 100 (see e.g. |20H ). Taking t ~ R c , this 



corresponds to wave numbers larger than L g = 27r(7 g G/i) . It is clear, from 
the physical interpretation of L* discussed in Section that gravitational 
damping can be neglected in the calculation of Vbh providing that L g > L*. 
This inequality gives 

Gfi < (27r 7 - 1 ) 5 / 3 (16vrV J B) 1 / 3 « (U^f^a'B) 1 / 3 . (Bl) 



Typically, the right hand side of ( |B1|) will be of order 10 -2 or larger, so this 
is not a very strong condition. 

Similarly we can consider the effects of friction. The dominant contribu- 
tion to friction comes from Aharonov-Bohm scattering of ambient particles 
off the string in the radiation dominated era (see Ref . |HJ ) . Assuming a situ- 
ation in which the wavelength of the perturbations is conformally stretched 
by the expansion of the universe (as it is in the present case), it is shown 
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in Ref. p2| that friction can only be important for wavelengths smaller than 
Xf = 7/(G/iT) _1 . Here T is the temperature and 7/ is a numerical coefficient 
of order one. The physical reason is that for A > A/, the friction term in 
the equations of motion will "switch off" before the perturbations cross the 
cosmological horizon and start to oscillate. 
The wave number corresponding to A j is 

L f = 2^] 1 RTG{i = ^ X GMT C , 

where T c is the temperature of the Universe at the time t c when the loop 
crosses the horizon, M is the mass of the loop at t c (which is also the mass of 
the resulting black hole) and we have used that the product RT is indepen- 
dent of time before the loop crosses the horizon. Using T c « (10G7V)~ 1/4 t- 1/2 
and t c ~ R c = M(2nfi)~ 1 (where M ~ 10 2 is the effective number of massless 
degrees of freedom) we have 



L, EE (lO^Vt 2 ^) 172 

\ m p / 



1/2 



Similarly to the case of gravitational radiation, friction can be neglected in 
the calculation of Vbh providing that Lf > L*. This condition is equivalent 
to 



Therefore, the range of values of G/i for which friction can be ignored depends 
on the mass of the black hole that one wishes to consider. 

Putting ( |BT| ) and ( |B2"D together and taking 7^ ~ l,j g ~ 10 2 , we conclude 



{a 2 B y2/9 ( ^\ 5/9 < Gfl < g . 10 -2 (a 2 s) l/3 



that damping is unimportant for values of G/i in the range 

M J 

Figure captions 

• Fig. 1 The function Vbh{L*) can be computed numerically from Eq. (|TT|) 
to arbitrary precision. We see that, essentially, the probability of black 
hole formation decays as an exponential function of L*. 
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Fig. 2 We represent the boundary of the region excluded by the con- 
straint (|22]) in parameter space (aGfj, versus B) for N = 27 (dashed 
line), N = 28 (solid line) and N = 30 (dotted line). The result is not 
very sensitive to the value of N and, quite generically, the bound is 
satisfied provided that aG/i~10~ 4 . 

Fig. 3 A plot of the number of loops N >e that would be seen in the 
sky at an angular size 0. 

Fig. 4 Schematic plot of the region excluded in parameter space by the 
observational constraints on loop nucleation. The horizontal boundary 
at Gfi ~ 10~ 6 is due to distortions in the CMB caused by loops with 
angular size larger than 10°. The slanted part of the boundary is due 
to 'unresolved' loops, with angular size less than 10°. The horizontal 
boundary at Gfi ~ 10~ 4 comes from limits on the abundance of pri- 
mordial black holes. The thick line is meant to represent values of the 



parameters that roughly satisfy the normalizations ( pop and (|36|), and 
which may be adequate for structure formation. 
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